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Abstract 

Due to the limited statistics so far accumulated in the Higgs boson search at the LHC, the 

D ■ Higgs boson property has not yet been tightly constrained and it is still allowed for the Higgs 

boson to decay invisibly to dark matter with a sizable branching ratio. In this work, we perform a 

^ \ comparative study for the Higgs decay to neutralino dark matter by considering three different low 

Q^ ■ energy SUSY models: the minimal super symmetric standard model (MSSM), the next-to-minimal 

^D ■ super symmetric standard models (NMSSM) and the nearly minimal super symmetric standard 

f^ \ model (nMSSM). Under current experimental constraints at 2(T level (including the muon g — 2 

(N ■ 

and the dark matter relic density), we scan over the parameter space of each model. Then in the 

allowed parameter space we calculate the branching ratio of the SM-like Higgs decay to neutralino 
^^ , dark matter and examine its observability at the LHC by considering three production channels: 

the weak boson fusion VV — )• h, the associated production with a Z-boson pp — )• hZ + Y or a pair 
of top quarks pp — )■ hti + X. We find that in the MSSM such a decay is far below the detectable 
level; while in both the NMSSM and nMSSM the decay branching ratio can be large enough to be 
observable at the LHC. We conclude that at the LHC the interplay of detecting such an invisible 
decay and the visible di-photon decay may allow for a discrimination of different SUSY models. 

PACS numbers: 14.80.Da,14.80.Ly,12.60.Jv 



I. INTRODUCTION 



As a cornerstone of the standard model (SM) and also the last undiscovered piece, the 
Higgs boson has been intensively searched in collider experiments. The foregone colliders 
LEP II and Tevatron yielded null search results, setting a lower bound of 114.4 GeV on 
the Higgs mass [1] and excluding a Higgs boson with a mass around 2Mw 2|], respectively. 
The ongoing Large Hadron Collider (LHC) took over the Higgs-hunting task and recently 
reported its search results. Based on an integrated luminosity of 4.9 fb~^ collected at 
^/s = 7TeV, the two experimental groups at the LHC independently further narrowed down 
the Higgs mass region (at 95% C.L. the CMS collaboration excluded 127-600 GeV while the 
ATLAS collaboration excluded 112.9-115.5 GeV, 131-238 GeV and 251-466 GeV) and both 
hinted to a Higgs boson around 125 GeV J3|. Such a finding has stimulated some theoretical 
studies for a Higgs boson near 125 GeV in low energy supersymmetry j^ and other models 



Of course, if the LHC hint of a 125 GeV Higgs from the di-photon channel is confirmed 
in the future, it would severely constrain or exclude those new physics models in which some 
new exotic decay modes (such as decaying invisibly into dark matter) are open and the di- 
photon rate is suppressed. But so far the statistics at the LHC is too small to confirm such a 
Higgs, let alone the precision measurement of the Higgs decay branching ratios. Therefore, 
experimentally it is still allowed for the Higgs boson to decay exotically, such as invisibly to 
dark matter, with a sizable branching ratio ^. 

Theoretically, the Higgs decay to dark matter can indeed occur in some new physics 
models, such as the gauge singlet extensions of the SM J7|, the SM with a heavy fourth 
generation |8|, the large extra dimension model [9|, the technicolor model [10'], the spon- 
ianeously broken R-parity models Jll| and the non-linearly realized supersymmetric model 
12| and the MSSM with a singlet [l3|. In this work, we perform a comparative study for the 
Higgs decay to neutralino dark matter in low energy SUSY by considering three different 
models: the minimal supersymmetric standard model (MSSM) 14J-ll6|. the next-to- minimal 
supersymmetric standard models (NMSSM) 17|,ll8l] and the nearly minimal supersymmetric 



"'^ In [6[ the authors used the hmited LHC statistics of the di-photon signal rates to set constraints on the 
invisible Higgs decay in the MSSM and found that the invisible Higgs decay branching ratio around 10% 
is allowed. 



standard model (nMSSM) [19|-l2l|. As will be shown, in both the NMSSM and nMSSM, the 
SM-like Higgs boson can decay to neutralino dark matter with a sizable branching ratio. 

In case that the Higgs boson decays to dark matter with a sizable branching ratio, de- 
tecting such a decay at the LHC will be important because in this case the conventional 
visible decays into 77, 66, tt^ WW'^*'^ and ZZ'^*'> are often suppressed. Obviously, the main 
production channel via gluon fusion gg ^ h is not usable because it just gives missing en- 
ergy. It was found through Monte Carlo simulations that the production via vector boson 
fusion ( VBF) pp — > hqq' and the associated productions pp — )■ hZ and pp — )■ hti can offer 



the opportunity to detect the Higgs decay to dark matter 22h26||. So in this work we choose 
these three production channels to display the observability of Higgs decay to dark matter 
in low energy SUSY. 

Note that although in the literature the Higgs decay to nutralino dark matter has been 
discussed in some specific model like MSSM, it is necessary to give a revisit in low energy 
SUSY for two reasons: (i) Different SUSY models usually give rather different phenomenol- 
ogy and it is interesting to perform a comparative study for different models; (ii) We want 
to know in the SUSY parameter space allowed by current experiments whether or not the 
Higgs decay to nutralino dark matter is detectable at the LHC. If in some model this decay 
is found to be accessible at the LHC, we further want to know how large the parameter 
space can be covered by searching for such an invisible decay at the LHC. 

This work is organized as follows. In Sec. II, we briefly describe the three supersymmetric 
models. In Sec. HI, through a scan over the parameter space we present the branching ratio 
of Higgs decay to neutralino dark matter and show the observability at the LHC. Finally, 
the conclusion is given in Sec. IV. 

II. THE SUSY MODELS 

In a renormalizable supersymmetric field theory, the interactions and masses of all par- 
ticles are determined by their gauge transformation properties and superpotential. The 
superpotential is a holomorphic function of chiral superfields $j D (0j,'?/'j, Fj), with (/),-, tjji 



and Fi being respectively the bosonic, fermionic and auxiliary fields, and takes a form 15| 



W = L%, + ^M*^*,!-, -f- ^y'^^^i^^^k (1) 



where the parameters Lj should be of dimension [mass]^ and is only allowed if ^j is a gauge 
singlet. The mass matrix M*-' can only be non-zero when the supermultiplets $» and $j are 
conjugates of each other under gauge transformation. And the massless coefficients y^^'' can 
only be non-zero when $j$j$jfc formed a gauge singlet. 

The MSSM is the most economized realization of supersymmetry in particle physics, 



which has two Higgs doublets H^ and Hd and its superpotential is given by [15| 

^MSSM ^^^^ ^fy^ . ^^^ (2) 

with Wp given by 

Wf = uYuQ ■ Hu - dYdQ ■ Hd - eY,L ■ Hd. (3) 

The MSSM has the so-called /i-problem, which can be solved in some extensions by intro- 
ducing a Higgs singlet 5*. Among these extensions the most popular ones are the NMSSM 



and nMSSM, whose superpotentials are [l7|, ll9| 



^NMSSM = ly^ + XSH^ ■ Hd + ^^3; (4) 

y^nMSSM = ^^^ xSH^ ■ Hd + ^pMlS . (5) 

The scalar potential in the Lagrangian contains the so-called F-term, D-term and soft-term 



m- 



where 



V'''^ = VF + Vn + Vsofu (6) 

5W 

Vp = F*'Fi, F*' = -W = ^; (7) 

Vn = ^i\Hd\' - \H^\r + 4iW\H^\' - \H^ ■ Hd^ (8) 

o z 

K^f ^"^ = ^hJ^«I' + ^hJ^-^I' + (B,H^ ■ Hd + h.c.)- (9) 

yNMSSM ^ rhl^\H^\'+ml^\Hd\''+ml\S\' + {XA^SH^ ■ Hd+^A^S' + h.c.) (10) 
yuMSSM ^ rhlJ\H^\' + ml^\Hd\^ + ml\S\^ + {XA^SH^ ■ Hd + isM'J + h.c.) (11) 

Here the parameter G is defined as G^ = gf + (^l with gi and g2 denoting respectively the 
U{1)y and SU{2)l couplings. With the Higgs fields H^, Hd and S developing respectively 



the VEV Vu Vd and f^, they can be rewritten as 
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V2 



(12) 



In both the NMSSM and nMSSM we have five complex scalar fields or ten real scalar degrees 
of freedom, whose mass eigenstates are obtained as 
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Here the three Goldstone bosons G^ and G^ will be eaten by the weak gauge bosons Z and 
W^ respectively. Then we have seven Higgs bosons, among which hi, /i2 and h^ are CP-even 
(with the convention rrih-^ < nih^ < i^hz), cii and 02 are CP-odd (ordered as m^^ < m^j), 
and H^ are the charged ones. 

In the NMSSM and nMSSM there are five neutralinos (x?)? which are the mixture of bino 
(B), wino (W^), higgsino {Hu, Hd) and singlino {S): 
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We assume the lightest neutralino is the lightest supersymmetric particle (LSP) and make 
up of the cosmic dark matter. 



For the purpose of our numerical analysis, we present the interactions o 
hi with bb, TT, WW and XiXv Iii the the NMSSM, they are given by 

rub 
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the Higgs bosons 

(15) 
(16) 

(17) 



(Sallfl + S,2nll + SrsUli) - V2kS,sNi,N: 



(18) 



where 11*/^ = NuNij + NijNu. The corresponding couphngs in the nMSSM can be obtained 
by setting k equal to zero (for MSSM, setting k, A, 5*13 and Ni^ to zero). 

For the NMSSM, in the basis x° = {.B,W'^,Hu,Hd,S), the tree-level neutralino mass 



matrix takes the form 17|, l28| 
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The corresponding mass matrix for the MSSM can be obtained by taking the upper 4x4 
matrix from the above expression, and for the nMSSM it can be obtained by set the term 
IkS to zero. When |/ie//| or Mg^ ^ M^, the lightest neutralino in the MSSM becomes 
bino-like, with a mass given by 



15| 
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(20) 



/i2 - Mf 
n the nMSSM the lightest neutralino is singlino-like and its mass can be approximated as 



20| 



A.1 
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(21) 



III. NUMERICAL RESULTS AND DISCUSSIONS 



We scan over the parameter space of each model under current experiment constraints, 
and for each survived sample we calculate the Higgs spectrum, decay branching ratios and 
production rates at the LHC. In our calculation we use the package NMSSMTools 



Hdecay 



and 



extend it to the nMSSM 2l|] . For the calculation of /i — t- 77 in the SM, we use the package 



29| . For the Higgs production cross sections, we use the code on the website [30| (this 



code is aimed at the MSSM, and we extend it to the NMSSM and the nMSSM). For parton 
distributions we use CTEQ6L [3l| with the renormalization scale and the factorization scale 
chosen to be the sum of the masses of the produced particles. 

In our scan we require the models to explain the cosmic dark matter relic density measured 
by WMAP [32] and also explain the muon anomalous magnetic moment at 2cr level. In 
addition, we consider the following experimental constraints: 



1. The LEP bounds on sparticle masses and on the Higgs sector from e~^e~ — )■ hZ{hA) 
followed by Z — 7- i~^i~ ,XiXiJ h ^^ hh and t^t^ [33]; We also consider the LEP- 
I constraints on the invisible Z decay ^, i.e., r(Z — t- x\xT) < ^-^^ MeV, and the 
LEP-11 constraints on neutralino production cr(e+e~ — > XiXi*) < 10^^ pb {i > 1) and 



a(e+e-^xM)<10"'pb(«>J>l) 



33|. 



2. The Tevatron bounds on sparticle masses and on stop or sbottom pair production 



followed by t — )• ex, blvi and b -^ bx (S^ |35| ; 



3. The recent LHC bounds on the Higgs sector from the measurement of the signal rates 

of 77, TT, WW^*^ and ZZ^*^ MM; 



4. The constraints from B-physics, such as 6 — )■ 37 and Bg — ?■ fifi 38| : 



5. The electroweak precision observables like Myy, sin^ 9^ and R^ 



In order to reduce the number of free parameters, for the gaugino masses we assume the 
grand unification relation 3Mi/5ai = M2/a2 = M^/a^ and thus we have only one gaugino 
mass parameter (we choose M2 in our calculation). In order to explain the muon anomalous 
magnetic moment at 2cr level for moderate tan/3 (< 20), for the smuon sector we assume 
the soft-breaking parameters to be 100 GeV 2l|. For other soft-breaking parameters in the 
squark and slepton sectors, we assume them to be 1 TeV except that for the MSSM we allow 
the third-generation squark mass parameters rrtg to vary in a wide range. The parameter 
Mn in the superpotential of the nMSSM is also fixed to be 1 TeV. Other parameters are 
scanned in the following ranges (among which Ma is the mass of the cos /?</)„ + sin Pipd field, 
i.e., the diagonal element of the doublet in the CP-odd Higgs mass matrix): 

1. For the MSSM: 1 < tan/3 < 20, 100 GeV < /x < 600 GeV, 10 GeV < M2 < 200 GeV, 
500 GeV < Ma < 3 TeV, 100 GeV < Mg < 2 TeV, -3 TeV < At^b < 3 TeV. In this 
space we scan two million random points and about three thousands points survived 
the experimental constraints. 



^ Note that the constraints from such an invisible Z decay are stringent for a wino-like or higgsino-Hke 
neutraUno, but become quite weak for a bino-Hke or singlino-Hke neutraUno. 



2. For the NMSSM: 0.1 < A < 0.7, 0.1 < /t < 0.5, 1 < tan/3 < 4, 100 GeV < (^u, M„) < 
1 TeV, 50 GeV < Ma < 150 GeV, < Ax < 1 TeV, -500 GeV < A^ < 0. In this 
space we scan ten milhon random points and about one thousand points survived. 

3. For the nMSSM: 0.1 < A < 0.7, 1 < tan/3 < 10, 50 GeV < (/i,M2,M«) < 1 TeV, 
< Ax <l TeV, < Mg < 500 GeV, -1 < ^i^ < 1. In this space we scan one billion 
random points and about 2 thousands points survived. 

Note that our scan is not a general scan over the entire parameter space. Since our purpose 
is to figure out if the invisible decay of the Higgs boson is accessible at the LHC, we only 
scanned over a part of parameter space which is potentially able to allow for a light neutralino 
and hence the Higgs can decay into the neutralino pair. 

In Figs. [1] and [2] we display the scatter plots of the parameter space which survive all 
constraints. Fig. [1] shows the Higgs decay branching ratios and the LHC (7 TeV) di-photon 
rate, while Fig. [2] shows the observability of the decay h — )■ XiXi through three production 
channels at the LHC (7 TeV). The 2cr sensitivity of the A' 



LAS detector shown in Fig. [2] is 
obtained from a Monte Carlo simulation carried out in |23l-l25| . For the decay h — )■ XiXi? the 
signature of the production via vector boson fusion VV -^ his two far forward and backward 
tagging jets of moderate pt with considerable missing momentum j/^ in the central region. 
For the productions pp — )■ hZ and pp -^ hti, the signatures are obvious: for the former it is 
two isolated high pt leptons from the Z-boson decay and large missing ^j. from the Higgs 
decay; for the latter it is di-leptons (or lepton plus jets) and large missing j/rp. 
From these figures we obtain the following findings: 

• In each model the SM-like Higgs can have a mass near 125 GeV, as hinted by the 
recent LHC results. 

• In the MSSM the SM-like Higgs boson dominantly decays to bb (just like in the SM), 
the decay h — )■ XiXi has a very small branching ratio (below about 10%), and the 
di-photon signal rate is close to the SM value. Due to the small branching ratio, the 
Higgs decay h -^ XiXi is far below the detectable level. 

• In the NMSSM the decay hi — )■ XiXi can be comparable to hi — )■ bb. In the region with 
a sizable decay ratio of hi — > XiXi; the lightest neutralino Xi is rather light (below 
hi/2) and the couphng of hi to Xi is large (see Eq. [18]). The diphoton signal rate can 
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FIG. 1: The scatter plots of the parameter space which survive all constraints listed in the text. In 
the upper frames, the samples denoted by crosses (sky-blue) show the branching ratio of /i — t- XiXi 
while the samples denoted by circle (magenta) show the branching ratio of /i — t- bb. The solid 
curves (blue) denote the SM prediction for the branching ratio of /i — )• bb. In the lower frames, the 
samples denoted by times (red) show the ratio (Jsusvipp -^ h ^ n)/<ysM{pp -^ h ^ 77) at the 
LHC (7 TeV). In this figure and below, ''hsM^ denotes the Higgs boson in the SM, 'h' denotes the 
lightest nutral Higgs boson in the MSSM, and '/ii' denotes the SM-like Higgs boson in the NMSSM 
and nMSSM (the doublet component of hi is over 



be sizably deviate from the SM prediction, either enhanced or suppressed significantly. 
In a large part of the parameter space, the Higgs decay hi — )■ XiXi is accessible at the 
LHC. 
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FIG. 2: Same as Fig.l, but showing the quantity ^^'^^'^ x Br{hi — ;• XiXi) with asusY {o'sm) being 
the SUSY (SM) Higgs production rates for the processes VV -^ h, pp ^ hZ and pp — )• htt. The 



sohd curves show the 2a sensitivity 



2a425|] of the ATLAS detector at the LHC (7 TeV) with 10 



fb ^,30/5 ^ and 100 fb ^ (the region above each curve is the observable region). 

• In a major part of the parameter space in the nMSSM, the decay hi — )■ x'JXi is 
dominant over hi — t- bb and thus observable at the LHC. The reason is the lightest 
neutralino is singlino-like and is always light, as can be seen from the neutralino mass 
matrix in Eq. (I2T1) . Also, from Eq. (ITSl) we see that the coupling ghix°x° '^^^ ^^ large 



10 



(near unity). Due to the new sizable decay hi — )■ x^Xi? the total width of the SM- 
like Higgs is greatly enlarged and thus its di-photon signal at the LHC is severely 
suppressed. So, if the recently observed di-photon signals at the LHC is verified in the 
near future, this model will be excluded. 

IV. CONCLUSION 

We examined the Higgs decay to neutralino dark matter in low energy SUSY by con- 
sidering three different models: the MSSM, NMSSM and nMSSM. We considered current 
experimental constraints at 2a level (including the muon g ~ 2 and the dark matter relic 
density) and scanned over the parameter space of each model. Then in the allowed param- 
eter space we calculated the branching ratio of the SM-like Higgs decay to neutralino dark 
matter and examined its observability at the LHC by considering three production channels: 
the weak boson fusion VV — )■ h, the associated production with a Z-boson pp — )• hZ + X or 
a pair of top quarks pp — )• hti + X. Our findings are: (i) In the MSSM such a decay is far 
below the detectable level; (ii) In the NMSSM it is accessible in a sizable part of parameter 
space; (iii) in the nMSSM it is detectable in a major part of the parameter space, (iv) When 
this invisible decay is sizable, the visible di-photon decay is suppressed. So, we conclude 
that at the LHC the interplay of detecting such an invisible decay and the visible di-photon 
decay may allow for a discrimination of different SUSY models. 
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